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TBE DYNAMIC-RESPONSE CHARACTERISTICS OF A 35° SWEPT- WING AIRPLANE AS DETER- 
MINED FROM FLIGHT MEASUREMENTS ^ 

By William C. Tetplbtt, Sthabt C. Bhown, and G. Allan Smith 


SUMMARY 

The longUvdiTwl and lateraJ^irectioTMl dynamie-respome 
cJuiracteristics of a 35° swept-wing JigMer-type airplane 
determined from flight measurements are presented and compared 
with predictions based on theoretical studies and wind-tunnel 
data. Flights were made at an altitude of 35,000 feet covering 
the Mach number range of 0.50 to 1.0 4 .. A limited amount 
of lateral-directional data were also obtained at 10,000 feet. 
Thefliglds consisted essentially of recording transient responses 
to pilot-applied pulsed motions of each of the three prirriary 
control surfaces. These transient data were converted into 
freguency-response form by means of the Fourier transformation 
and compared wilh predicted responses calculated from the 
basic equations of motion. The equations, or tranter func- 
tions, that best describe the various measured responses were 
evaluated by a curve-fUting process involving the use of tem- 
plates and an analog computer. By this method it was generally 
possible to find equations, of simple form, that closely matched 
the experimental frequency responses between 1 apd 10 radians 
per second and at the same time adequately described the re- 
corded time histories. 

Experimentally determined transfer functions were used for 
the evaluation of the stability derivatives that have the greatest 
effect on the dynamic response of the airplane. The values of 
these derivatives, in most cases, agreed favorably with predictions 
over the Mach number range of the test. 

INTRODUCTION 

In the design of automatic-control equipment for high 
performance aircraft, the dynamic response characteristics 
of the aircraft must he considered. It is desirable to express 
these characteristics as transfer functions which are expres- 
sions that describe the motion of the airplane for the various 
flight conditions of interest. The airplane can then he rep- 
resented as a single element in a more complex closed-loop 
system. 

Often these dynamic characteristics can be predicted by 
using stability derivatives obtained from wind-tunnel tests. 
In many cases, however, particularly in the transonic speed 
range, flight-test procedures are desirable to document the 
dynamic behavior of the airplane. Flight teste also serve 
the additional purpose of enabling comparisons to be made 
with predicted results, thus aiding in the development of 
more reflned prediction methods. 


This report describes the results of a flight investigation 
in which the dynamic-response characteristics of a 35“ 
swept-wing airplane have been evaluated through the Mach 
nmnber range of 0.50 to 1.04. Responses to transient 
rather than sinusoidal control inputs have been chosen for 
analysis because of convenience in making flight measure- 
ments. Certain useful information can be obtained directly 
from data recorded in transient-tune-history form; however, 
when converted to frequency response form, the dynamic 
response of the aircraft is presented in a standard manner 
independent of the particular input used, and the dynamic 
characteristics are more readily apparent. The frequency- 
response data are also applicable to analysis by conventional 
servomechanism methods. ; 

Experimental frequency responses be described by 
analytical expressions termed “transfer functions” which are 
evaluated by a method utilizing a set of frequency-response 
templates and an analog computer. Whenever possible, the 
coefficients of these transfer functions are expressed in 
conventional stability-derivative form and compared with 
wind-tunnel and theoretical estimates. ' 

NOTATION 
lift coefficient 
rolling-moment coefficient 
pitching-moment coefficient 
yawing-moment coefficient 
side-force coefficient 

the operator, ^ 

moments of inertia about the X, Y, and Z axes, 
slug-ft* 

product of inertia, slug-ft^ 

Mach number 

real and imaginary parte of a complex quantity 
wing area, sq ft 
velocity, ft/sec 
weight of airplane, lb 

wing span, ft ■ . 

mean aerodynamic chord, ft 
acceleration due to gravity, ft/sec*. ■ 

V=1 

mass of airplane, slugs 

normal acceleration, V(J)L—g), ft/se(^ (except as 
noted) 
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Sapersedes recently declassified NAOA. RM A51Q37 by WlHiam 0. Triplett and G. Allan Smith, and NAOA BM A52U7 by w nuam O. Triplett and Stuart 0. Brown* 
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rolling velocity, radians/sec 
pitching velocity, radians/sec 

radians/sec* 

dynamic pressure, Ib/sq ft 
yawing velocity, radians/sec 
time, sec 

angle of attack, radians (except as noted) 
radians/sec 

sideslip angle, radians (except as noted) 
flight-path angle, deg 

control deflection, radians (except as noted) 
total aileron deflection, radians (except as noted) 
elevator deflection, radians (except as noted) 
rudder deflection, radians (except as noted) 
damping ratio 

root of the characteristic equation 

real part of a complex root 

angle of bank, radians 

phase angle, deg 

angle of yaw, radians 

frequency, radians/sec 

natural frequency of oscfllation, radians/sec 

imdamped natmtd frequency, radians/sec 

S(^T^ per radian 
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^0, 
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op 
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TEST EQUIPMENT 

The test airplane was a standard North American F-86A-6 
airplane with external instrument booms added as shoivn in 
figure 1. The physical characteristics of this airplane are 
described in table I. 
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TEST PBOCEDUEES 

The flight-test procedures consisted of recording airplane 
responses to pulse-type disturbances of each of the three 
control surfaces separately. Flights -were made at an altitude 
of 35,000 feet through liie Mach number range of 0.50 to 
1.04. The trim lift coefficient varied from 0.51 to 0.12. A 
limited amount of data were also obtained at an altitude of 
10,000 feet over this range. 

During each test run one control surface was deflected and 
then returned to the trim position to provide a pulse input 
while the other two surfaces were held flxed. After applica- 
tion of the pulse, all three surfaces were held flxed until the 
oscillatory motion of the aircraft had essentially subsided. 
Sample time histories of airplane responses to the three 
control-surface inputs are shown in figures 2 and 3. The 
variation with Mach number of trim angle of attack and 
trim elevator angle are plotted in figure 4. 

AU flight runs at speeds below a Mach number of 0.95 
were initiated from trimmed level flight, but to obtain data 
at the higher speeds it was necessary to dive the airplane. 
Altitude changes up to 2,000 feet were encountered during 
each diving run although there was little variation in Mach 
number. For analysis purposes the altitude and dynamic 
pressure were assumed constant at their average values 
during each run. 



Standard NACA instruments were used in measuring the 
following quantities: Pitching, rolling, and yawing velocities 
were measured by rate gyros with direct optical recording. 
Normal acceleration was measured by an air-damped vane- 
type accelerometer. Elevator, rudder, and aileron deflec- 
tions were measured by control position transmitters that 
were linked directly to the control surfaces, and were re- 
corded on an oscillograph. An^es of attack and sideslip 
were measured by a vane-type pickup and recorded on an 
oscillograph. True Mach number and altitude were obtained 
from the nose-boom airspeed system described in reference 1. 
All recordings were synchronized at 0.1-second intervals by 
a common timing circuit. 

The rate gyros used to measure angular velocities and the 
accelerometer each had a sufficiently high natural frequency 
so that corrections to the data for instrmnent dynamics 
were not required. Tests of the control position recorders 
also indicated negligible dynamic lag through the frequency 
range of interest. 




FiauBE 2, — Sample flight records of normal acceleration, pitching 
velocity, and elevator angle at a Mach number of 0.81. 
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(a) Rudder input. 

Figtjbb 3. — Sample flight records of yawing velocity, roUing velocity, 
and sideslip angle at a Mach number of 0.81. 


METHOD OF ANALYSIS 

The procedures used in determining the airplane response 
-characteristics from the measured transient time histories 
•can be divided into three distinct steps which are outlined 
in the following sections. These are: frequency response 
•calculations, evaluation of transfer functions, and a final 
•check utilizing an analog computer. Wherever possible, 
^stability derivatives are extracted from the coefficients of 
the transfer fimctions. 


FBEQDENCT EESPONSE CALCULATIONS 

The first step involves the conversion of transient time 
histories into the frequency domain. 

Under certain conditions a tune function such as the 
pitching velocity g{t) can he transformed into a complex 
'frequency function 2(^) ky means of the Fourier integral 
a*elation 






This integral must be evaluated from time zero to infinity 
for each frequency a at which g^{iu) is desired. Obviously, 
the Integration can be accomplished only if the behavior of 
g(f) is known for an infinite time. Since g(f) can be measured 
only for a finite time T, it is necessary for the system to 
reach steady-state conditions before time T such that g(0 
may be expressed analytically between the time limits T 
and infinity. In addition, the product must con- 

verge as t approaches infinity. 

Examination of the data used in this analysis showed that 
in every case the responses and the forcing fimction became 
constant after a finite time interval T. In the case of the 
pitching-velocity response to an elevator input these steady- 
state values may be expressed as g(T) and 5d(T), respectively, 
and thus equation (1) may be divided into two parts as 
follows: 

g('Uo)=J' gff)e~*“‘dt+g^T)J 
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(b) Aileron input. 

Figubb 3; — Concluded. 
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FiaxmB 4. — ^Variation of trim angte-of-attack and trim elevator angle 
with Mach number. 


Tho second integral can be evaluated analytically so that 


qiiw)=& r 

W Jo 


( 2 ) 


For computational purposes equation (2) may be divided 
into its real and imaginary parts represented by the s5Tnbols 
R and I, respectively, so that 


where 




R=-^sm 
7=-^ cos 


cos o)t dt 
sin cot (ft 


Tho two integrals involving the transient part of q{() may 
may be evaluated by any of several approximation methods. 
In analyzing the data of this report, Simpson’s rule was used 
to find the area under the product curves g(t) cos cot and 
g(t) sin cot. Values of g(t) were tabulated at 0.05-second 
intervals and tho integrations were carried out at each 
integral value of frequency from 1 to 10 radians per second. 

After obtaining the numerical values of R and I at each 
frequency, g(fco) can be expressed in polar notation such that — 

g(ico) = |g|e**« 


where the omphtude, |2|= angle, 
arc tan IjR. 

Tho above integration process was repeated for the elevator 
forcing function Se(t) to determine 


The ratio of gJS, is then expressed as 



where |g/5„l is the amplitude ratio and is the phase 
difference between the two quantities, 

Similar calculations were made for the normal acceleration 
response to an elevator input and for the rolling, yawing, and 
sideslip responses to both rudder and aileron inputs. 

In order that flight-test results can be correlated directly 
with wind-tunnel measurements it is sometimes desirable to 
correct the angular velocities measured in flight so that they 
conform to the system of stability axes rather than the body 
axes about which the recording instruments are alined. The 
necessary correction is described in Appendix A. 

A further consideration in frequency-response calculations 
is the shape of the forcing fimction. As discussed in Appendix 
B, the type of input used in flight definitely places a limit on 
the accuracy of the Fourier analysis. In general, to obtain 
the widest usable frequency range, a pidse input shoifld be 
used. When low frequencies are desired, a step input is 
preferable, although this type of disturbance may result in 
motions that exceed the ranges of linearity. 

The final resiilt of the above calculations is a graphical 
representation of the airplane transfer function (frequency 
response) for a given test condition. This is plotted as 
curves of amplitude ratio and phase angle versus frequency. 
The remainder of this section of the report is concerned with 
the determination of analytical expressions for these curves. 

Reference 2 is one of many sources that explain in more 
detail the use of the Fourier transform for problems of this 
type. 

DETBEMINATION OF TBANSFEE-FDNCmON COEFFICIENTS 

Dynamic-response templates. — ^A graphical method was 
used to find the type of transfer function that beat defines 
each of the measured frequency responses. Numerical 
values of the transfer coefficients were also determined in 
the same operation. This method involved the use of a set 
of templates developed by C. S. Draper of the Instrumenta- 
tion Laboratory of the Massachusetts Institute of Tech- 
nology. These templates are described fully in reference 3, 
and their design is based on the following considerations: 

1. A rational function of complex frequency D can, in 
general, be expressed as 


F{D)=K 


{D+a^{D+a^ 

(P+h;){D+h^ 


(D+ga) 

iD+h„) 


(3) 


where the a* and ht may be either real or complex, and, if 
complex, always appear as conjugate pairs and where K is 
any real number. When is a real number, a factor 
(Z>-4-gi) may be expressed as 


i(l-t-rD) 

T 

where r=— • When and (i(+i are a complex conjugate pair 

di 

the two factors (D-l-(i,)(D-f can be expanded, by intro- 
ducing new coefficients, into the form (D^-l-2fco„D-l-co„^) 


418072— G7 77 





1222 


EEPOET 1260 — NAnONAL ADVISOBT COMMITTEE FOB ABBONACTICS 


where a^=afi,t+i and 2fwo=Ot+a<+i. J’actoring out 
then gives 

»»’(l+2f£+S) 

The terms — and that appear outside the parentheses, 

being real numbers, may be grouped together with the 
multiplying factor K. When an «« or 6* is zero then a factor 
D win be present in either the numerator or the denominator. 
Then, F(D) in factored form will consist entirely of combi- 
nations of the three forms. 

Since we are interested only in the steady-state frequency 
response, the complex variable D can be replaced by the 
frequency variable ia, (where i=-yj — 1) and each of the 
three factors may be written in polar form as follows: 

For 

where 

$i=tan~^ T<a 


( 4 ) 


for (l4-2f--f 
V Wi. W* / 




where 


and for D 
where 


!-(«*/«,*) 

uo=u>e**3 

$8==tian“' 00=90° 


(6) 


2. When tw is taken as a nondimensional frequency 
variable, all first-order terms of the form (1-f •n&j) can be 
defined by a single amplitude curve and a si n g le phase-angle 
curve as shown in figure 5. In a similar manner all possible 

second-order terms of the form l-f-2f — can be 

\«»/ 

defined by a family of curves with as the nondimensional 
frequency variable. There will be a different pair of curves 
for each value of damping ratio f, as shown in figure 6. For 
the factor ico (eq. (6)) the amplitude is simply equal to « and 
the phase an^e is a constant 90°. 

3. Since a typical transfer fimction may consist of inore 

than one of these factors in both numerator and denominator, 
it is of great advantage to plot the amplitudes to a logarithmic 
scale so that multiplication or division of factoi-s of the type 
shown in equations (4) to (6) may be accomplished by mere 
graphical addition of the amplitudes. The nondimensional 
frequency has also been plotted to a logarithmic scale as 
shown in figures 5 and 6. Then as frequency approaches 
either zero or infinity the amphtudes become asymptotic to 
straight lines. In both figures as frequency decreases the 
amplitudes approach asymptotically the value of unity. At 
high frequencies the amplitude of (l-brio) is essentially 
equal to to and thus it plots as a striaght line of slope 1 (fig. 
5). Similarly, the amphtude of a second-order term ap- 
proaches the value frequency increases and is 

asymptotic to a straight line of slope 2. Since a second-order 



TW 


FiauBB 6. — ^Variatioii of the first-order term (l+rtw) with^fnendl- 
menaonal frequency ra. 


term of this type normally appears in the denominator, it has 
been plotted as 

t+^4:+OT‘ 


and thus the slope of the amplitude curve is —2 (fig. 6). In 
both plots the asymptotes intersect at unity on the non- 
dimensional frequency scale. This intersection is termed 
“the breakpoint.” 

The phase angles are plotted to a linear scale because when 
factors of the types shown in equations (4) to (6) are com- 
bined, the resultant phase angle of F(D) is merely the alge- 
braic sum of the individual angles. It can be seen in figures 
6 and 6 that the angles of these first- and second-order terms 
approach 90° and —180°, respectively, at high frequencies. 
ThTis an amplitude slope of 1 corresponds to an angle of 90°, 
while a slope of —2 corresponds to an angle of —180°. In 
general, as w approaches infinity the amplitude of F(D) 
approaches where n and m are defined in equation 

(3). Then the amplitude slope is equal to (n — m) on the 
logarithmic plot and the phase angle is 90 (ti— m)°. 

The templates used in the analysis of the data were accu- 
rate representations of the curves shown in figures 6 and 6 
cut from tratoEparent material with the breakpoints marked. 
Twelve pairs of second-order templates were included in the 
set to give values of damping ratio from 0.1 to 1.0 in 
increments of 0.1 and in addition values of 0.06 and 0.16. 
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(a) Amplitude. 


Fiotjrb 6. — Variation of the second-order term 




with nondimensional frequency for various values of damping 
ratio. 


The computed frequency-response data (amplitude ratios 
and phase angles) were plotted to the same logarithmic scale 
as the templates, and then by a trial-and-error approach the 
template or combination of templates that best matched 
the given data was determined. To determine the proper 
combination, various amplitude templates were positioned 
in such a way that their algebraic sum matched the experi- 
mental ampHtude ratio curve. For each amphtude tem- 
plate the corresponding phase-angle template was placed 
on the experimental phase-angle curve so that the break- 
point of each pair was alined with respect to frequency. 
The various pairs of templates were then adjusted imtil the 
sum of the individual amphtudes and phase angles most 
nearly matched the test data. The frequency at which the 
brealq)oint of each pair occurred (wa) was noted and then 
the appropriate values of t and w, were determined directly, 
since for each first-order term 

_ 1 
(^B 

and for each second-order term 


0>rt — 

In general, each factor will have a different breakpoint. 

At first glance, this procedure may seem extremely tedious 


but it has been found that with practice one can determine 
coeflBcients quite rapidly in this manner. For a more com- 
plete description of the principles involved in this type of 
graphical representation reference 4 is recommended. 

Analog Computer. — ^As a final step in the calculations a 
Reeves analog computer was used to refine the previously 
determined values of transfer-fimction coefficients. This 
was accomplished by placing the actual time histories of 
the control motion into the computer by means of an input 
table. Then this input was fed into a circuit representing 
the equations of motion as obtained in the previous step and 
the outputs of the machine were obtained. These outputs 
were compared to the actual time histories of the airplane 
responses as obtained in flight. By changing dial settings 
on the computer the transfer coefficients could be adjusted 
until the output of the computer most nearly matched the 
actual flight data. (Such a comparison is shown in figs. 2 
and,o.)- Thus, in addition to providing a refinement of the 
coefficients obtained with the templates, this step also 
resulted in a check of aU previous calculations. In general, 
operations of this type can be conveniently accomplished 
on the computer only when the form of the transfer function 
is known. 



Fiqubb 6. — Concluded. 
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RESULTS AND DISCUSSION 

In the following paragraphs the longitudinal character- 
istics and the lateral-directional characteristics are discussed 
separately. Flight-evaluated frequency responses, transfer 
functions, and stability derivatives are presented and 
compared with predictions based on wind-tunnel data and 
theoretical studies. 

LONGITUDINAI. EESPONSE CHABACTEEIS'nCS 

Frequency response. — ^Plotted in figure 7 are typical 
flight-evaluated frequency responses of pitching velocity 
and normal acceleration to elevator inputs. The purpose 
of these figures is to show general trends with varying Mach 
number, and therefore smooth cmwes have been faired 
through the calculated test points. In most cases more than 
one flight record was analyzed at each flight speed and only 
at the highest speeds was there appreciable inconsistency. 
Each transient record was analyzed through the frequency 
range of 1 to 10 radians per second and some scatter in the 
data was noted at either end of the frequency range. For 
this reason, as well as for clarity, portions of some of the 
responses have been omitted from figure 7. 

The amplitude-ratio plots show the increase in natural 
frequency and the decrease in amplitude which reflect the 
changes in static margin and elevator effectiveness with 
Mach number. 

To afford a comparison with estimates made from wind- 
tunnel data, typical frequency responses for a Mach number 
of 0.81 have been replotted in figure 8. Shown here are the 
actual data points obtaiaed from a single transient input. 
Also shown as solid lines are responses predicted by use of 
the stability derivatives of table II (ref. 5) and transfer 
functions developed in Appendix C. The scatter between 
individual data points is typical for responses in the lower 
speed range; at Mach numbers near 1.0, however, there was 
somewhat more scatter in the individual data points. 

Transfer functions and stability derivatives. — ^The airplane 
transfer functions were determined from the frequency re- 
sponses of figure 7 by using the templates described in the 
Analysis Section. Final values of the coefficients were 
obtained from the analog computer by matching the original 
time histories. It was found that the airplane responses 
could be adequately defined by transfer fimctions of the 
form: 

ff OiJD-\-Got 

S,~D'+iD+k 

and 

n _ Go, 

srD^+hB+k 

As indicated in figure 2 the output of the analog computer 
closely matched the transient motions measured in flight. 
These equations are also of the same form as the theoretical 
transfer functions developed in Appendix C. 

The variation with Mach number of the coefficients 
and f are plotted in figure 9. The dampiog-ratio curve was 
faired in accordance with additional data which more clearly 
defined the sharp variations between Mach number of 0.88 


and 0.96. An indication of the consistency of the data is 
given^by the scatter of the pomts at each flight speed. 




(a) Pitching velocity response to elevator input. 

FiatnaB_7. — Longitudinal frequency responses at various flight Maoh 
numbers at an altitude of 36,000 feet. 




Frequency, w, radians/sec 

(b) Normal acceleration response to elevator Input. 
FiauEB 7. — Concluded. 




Phase ongle. 
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(a) Pitching velocity response. 

Fioueb 8. — Comparison of experimental and predicted longitudinal 
frequency responses at a Mach number of 0.81 (altitude 36,000 
feet). 




(b) Normal acceleration response. 
Fiqurb 8. — Concluded. 
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Figtoe 9. — Variation of longitudinal undamped natural frequency 
and damping ratio •with Mach number. 


Plotted in figures 10 and 11 are the stability derivatives 
Omai ^nd Thcse -nrere evaluated from the 

transfer coefficients as indicated ia Appendix C. Because of 
the difficulty in finding reliable values of the coefficient 
wind-timnel values of Cl^ (also shown in fig. 10) were 
used to calculate Comparisons are made with 

estimates sho'wn in table II and in the case of the results 
of this investigation are compared to unpublished static 
flight data. 

LATEBAL DIEECTIONAL BESFONSE CHABACTEBISTICS 

Freqtnency responses. — Plotted in figure 12 are typical 
flight-evaluated frequency responses of rolling velocity, 
yawing velocity, and sideshp angle to rudder and also to 
aileron inputs. These were all obtained at an altitude of 
35,000 feet; responses at 10,000 feet showed similar character- 
istics and have not been plotted. As in the case of the 
longitudinal frequency re^onses smooth curves have been 
faired through the calculated test points. 

The responses to rudder inputs have not been plotted at 
frequencies greater than 8 radians per second because there 
was considerable scatter and also a lack of well-defined trends 
in the data at the higher frequencies. The aileron responses, 
however, are shown to 16 radians per second. Wherever 
necessary for clarity or because of erratic data, parts of some 
of the curves have been omitted. The /3/5a response is shown 
at only three speeds because of a failure in the sideshp-angle 
recording system. 



Figubb 10. — Variation 'with Maoh number of 'the stability derivatives 


and h wit Mach number. 




Fig'obb 11. — ^Variation ■with Mach number of the stability deriva- 
tives Cx^ and with Mach number. 
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(a) RolUng velocity resjKDiise to rudder input. 

FiaUEB 12. — Lateral-directional frequency responses at various flight 
Mach numbers (altitude 36,000 feet). 


With minor exceptions, the curves show consistent and 
gradual variations with Mach niunber. One such exception 
can he seen in figures 12 (d) and 12 (e) where the amplitudes 
of pjSa and r/5a at a Mach number of 0.61 lack the customary 
resonant peaks. This is the result of time histories in which 
there was no oscillatory motion. This unusual characteristic 
can be explained by reference to the predicted transfer func- 
tions developed in Appendix C and discussed in the following 
section. The predicted p/5a response equation for a Mach 
number of 0.6 at 36,000 feet is 

V_ 21.1Z)(D>4-0.455l>-h6.74) 

5a (Z>-h0.00113)(Z>-|-2.203)(I>*-f-0.438I>-+-7.25) 

It can be seen that the two quadratic terms are nearly identi- 
cal and thus the oscillatory mode is effectively canceled. 



(b) Yawing velocity response to rudder input. 
Fiqube 12. — Continued. 


Since rolling and yawing motions are coupled, the r/5a re- 
sponse must exhibit the same characteristics at this particu- 
lar speed. 

Another interesting point with regard to figure 12 (e) is 
the wide variation in phase angles at different flight speeds. 
Predicted transfer fimctions indicate that at low speeds 
(below a Mach number of 0.7) where 0^^ is negative, the 

phase angles approach —270° asymptotically with increasing 
frequency. Unpublished windrtimnel data indicate that 
near a Mach niunber of 0.7 there is a transition in which 
(7,,^ becomes positive and consequently three of the co- 

efl&cients in the numerator of the transfer function change 
sign. The result is an increase of 180° in the high-frequency 
phase lag. 
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Frequency, <u, rod ians/sec 
(c) Sideslip angle response to rudder input. 
Figubb 12. — Continued, 



(d) Bolling velocity response to aileron input, 
Figtoe 12. — Continued. 
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The frequency-response test points derived from flight data 
at a Mach number of 0.81 have been replotted in figure 13 
which shows all six responses for the 35,000-foot altitude. 
These results are typical in indicating the degree of scatter 
usually encountered in the Fourier analysis of a particular 
flight record. Plotted as sohd lines for comparison are pre- 
dicted responses that have been calculated using estimates 
of the various stability derivatives presented in table II, 
which were obtained from reference 6 and also from wind- 
tunnel tests by the manufacturer. These calculations were 
made as shown in Appendix C, using the exact linear fourth- 
order response equations. The agreement between measiued 
and predicted responses is generally good ^cept, in some 
cases, at the extremes of the test frequency range. This is 




(a) Rolling velocity response to rudder input. 

Fiqube 13. — Comparison of experimental and predicted lateral 
directional frequency responses at a Mach number of 0.81 (altitude 
36,000 feet). 


particularly noticeable in figure 13 (d) where the measured 
frequency response of pjSa indicates a mode of motion not 
consistent with the rigid-body equations of Appendix 0. 
The increase in phase angle and the “up-tum” in the ampli- 
tude curve at high frequencies indicate the presence of an 
additional mode which is probably related to the primary 
bending frequency of the wing (about 8 cps). 



Frequency, -w, radions/sec 

(b) Tawing velocity response to rudder input. 
Fiqtjbb 13. — Continued, 
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Frequency, cu, rodions/sec / 


(o) Sideslip angle response to rudder input. 

Fioubb 13. — Continued. 

In addition to the pulse-type inputs, a step disturbance of 
the rudder was used in one instance at a Mach number of 
0.81 in order to define more clearly the low-frequency por- 
tions of the p/5r and r/3, frequency responses. Results of 
this analysis from 0.1 to 1.0 radians per second are plotted in 
figures 13 (a) and 13 (b) and verify the prediction of a sharp 
attenuation in amplitude of r/S, at a frequency of 0.4 radian 
per second. 

Also shown as dotted lines in figures 13 (b), (c), and (d) are 




Frequency, w, rodions/sec 

(d) Rolling velocity response to afleron input. 
Fiotteb 13. — Continued. 


responses computed from predicted transfer functions that 
have been simplified as indicated in the following paragraphs. 

THEORETICAL TRANSFER FUNCTIONS 

It is shown in Appendix C that the characteristic equation 
A can be factored into the form 

A=P(D-Xi) (P-Xs) (IP+CiD+Ca) 

where Xj and Xa are the spiral and rolling roots, respectively, 
and where Ci and Cs are coejB&cients that define the oscillatory 
mode. By neglecting \ (which is usually very small) and by 
omitting other small terms that appear in the numerators of 
the various response equations, three of the six responses may 
be reduced to the following simple forms: 

r_ BJ) 
dr D^-\-cJ)-\-Ci 

^ — Bz 

Sr 

V_ Az 
Sa D Xj 
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Frequency. <u. radions/sec 


(e) Tawing velocity response to aileron input. 
Figtjbb 13. — Continued. 




(f) Sideslip angle response to aileron input. 
Figtob 13. — Concluded. 


It is also shoAvn that by making additional assumptions as 
to relative magnitudes, the coefficients Xj, Ci, and Ct can be 
expressed as Lp, —(Nr+Yfi), and N/, respectively. Fur- 
thermore, since Bt^Ns/ and As^Ls^', equations (7) can be 
Tpritten as 


r 

Sr 


n>;d 

'W=W^WWW 


(8a) 


p -Ns/ 

5r D^-(Nr+T,)D+N/ 


(8b) 


P_ 

8a D — Lp 


(8c) 


It can be seen that in equations (8a) and (8b) the spiral 
and rolling modes are completely neglected, and yet, as shown 


in figures 13 (b) and (c), these simplified transfer fimctions 
yield responses that are almost identical to those obtained 
from the “exact” equations for frequencies greater than 1 
radian per second. Similarly, the response computed from 
equation (8c) closely matches the exact response (fig. 13 (d)) 
over the frequency range shown except that it omits the small 
peak normally associated with the oscillatory mode. The 
spiral mode which has been neglected in all three simple 
equations appears to have no effect on the calculated air- 
plane response except at frequencies well below 0.1 radian 
per second. 

EXPERIMENTALLY DETERMINED TRANSFER FUNCTIONS 

In the analysis of the fiight data it was found that the 
frequency responses of r/5r, iS/S,, and p/6a could be successfully 
simulated by simple transfer fimctions of the same forms as 
equations (7). Solutions of these equations on the analog 
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computer, using final “best” values of the numerical co- 
efficients with actual control motions as recorded in flight, 
resulted in outputs that closely matched the measured time 
histories of r, /3, and p as shown in figure 3. This fact implies 
that the modes of motion that are neglected in each case have 
very little effect on the time response to a pulse-type input. 

By use of the coefficients that best describe the measured 
time histories, frequency responses were calculated for com- 
parison with those derived directly from flight data. Ex- 
amples of these calculations are shown by the dash-dot lines 
in figures 13 (b) and (d). These curves, in general, match 
the experimental points closely for frequencies between 1 
and 10 radians per second. 

It is apparent from equations (8a) and (8b) that the same 
information can be obtained from either r/5, or /3/5,. Coef- 
ficients evaluated from each of these responses agreed favor- 
ably in most cases; however, because of indications that the 
yaw rate gyro possessed dynamic characteristics superior to 
those of the sideslip vane, only the yawing velocity responses 
were used in the final calculations. Natural frequency and 
damping ratio of this mode are plotted in figure 14. 




Fiourb 14. — Variation of directional undamped natural frequency 
and damping ratio with Mach number. 


The transfer functions of the three remaining responses 
p/5,, r/5a, and /3/5a were not amenable to simplification. 
However, it was found that the p/5, response could be 


matched satisfactorily by a transfer function of the type 

p (h{D+a^ (H+Os) 

5,~(H-Xs) (Z?*-i-CiZH-c) 


which is the same form as developed in Appendix C except 
that the spiral mode has been neglected. As written here, 
fli is identical to Li' while Oj and are complicated combina- 
tions of derivatives that cannot be readily simplified. Al- 
though this equation closely describes the measured time 
histories (fig. 3), it was difficult to find unique values of the 
numerator coefficients. Changes in one of these could be 
compensated for by corresponding changes in the other two, 
and the values were not considered to be reliable enough for 
presentation. 

Definition of the r/5a and /S/5a responses required fourth- 
order transfer functions that include all three modes and, 
because of practical difficulties involved, no attempt was 
made to evaluate the coefficients of these responses. 

STABIUrr DERTVATIVES 


In addition to the quantities L^, N,-\-Tp, N^', Ni/, and 
LtJ, that were determined as mentioned in the preceding 
paragraphs, the coefficient Lp was evaluated from the time 
histories of rolling and yawing velocity as outlined in Ap- 
pendix D. 

The quantities JV/s, L^, Ni^, and Ls^ can then be calculated 
from N/, Lfi, Ni/, and LtJ by usiog the foUowing^ expres- 
sions obtained from the relationships developed in Appendix 


C: 






U-rxN^' 

l—TxTz 


Ni 


r 


Ni/-rzLi/ 

l—TxTz 


L^a- 


LiJ-rxNiJ 

l—rxTx 


Because rx and are very small quantities, wind-tunnel 
estimates of ij/ were assumed to be sufficiently accurate 
to use in the calculation of iVj^. The term VxNiJ was 
completely neglected in evaluating Li^. 

Finally, from the definitions given in the notation it is 
possible to evaluate the derivatives (7,^, Gi^, 0„^ , and 



The anal3rBis methods used herein do not allow the separa- 
tion of the damping term Nr+Yf. As compared to N„ the 
term is small and can generally be predicted accurately 
from wind-tunnel measurements. Therefore values of Orp 
given in table II were used iu calculating fr6m the 


quantity Y^. 
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The flight evaluated derivatives for both altitudes are 
plotted against Mach number in figures 15 and 16. These 
are compared to the predicted values listed in table II. 
Tirfough the speed range of the test the predictions for both 
altitudes are essentially the same, except as noted in the 
plot of Ci^. 

The correlation between predicted derivatives and those 
evaluated from flight at 35,000 feet is generally good except, 
in some cases, at speeds near a Macii number of 1.0 where 
the predictions are apt to be inaccurate. Unpredicted varia- 
tions with altitude are also apparent in the flight values of 
Utp, and Un, . At a Mach number of 0.8 the value of 

Gn, for 10,000 feet is less than one half the value at 35,000 
feet. Fhght data of reference 6 when expressed in this form 
show a similar trend. The plot of C7„j^ indicates the same 

tendency to a lesser degree, while the 10,000-foot value of 
(at iU=0.8) is some 50 percent higher than the value at 
35,000 feet. A possible explanation for these discrepancies 
is that since dynamic pressure varies with altitude, structmul 
deflections due to aerodyoamic and inertial loads may cause 
changes in the effective values of the derivatives. 

Values of Ci^ determined in the present investigation agree 
favorably with wind-t unn el results, while those reported in 
reference 7 (obtained from static flight tests of the same air- 
plane) are much smaller in magnitude. It appears, however, 
that-the results of reference 7 are subject to error becaiise of 
^e simplifying assumptions made. A more rigorous ap- 
proach would have resulted in laig:er values of this derivative. 

Ejcamination of figure 16 shows the control effectiveness 
derivatives G„g and Gig to have similar variations with 

increasing Mach nmnbOT, and in each case the measurecd 
values are generally smaller than predicted. Values of (7,^^ 

obtained in the present investigation agree closely with those 
presented in reference 8 which again were evaluated from 
flight measurements of the same airplane. 

In this iavestigation there was no evidence of nonlinear 
variations of rolling or yawing-moment coefficients with 
p, r, or p. This was concluded because (1) the period and 
damping of the oscillations following a control input were 
essentially constant in every case (no systematic variations 
with amplitude), and (2) the experimental time histories 
coifld be matched, in general, by differential equations with 
constant coefficients. 

No conclusions are drawn as to nonlinear moment coeffi- 
cient variations with da or because the magnitudes of the 
control inputs were not varied appreciably during the. tests. 
ThSy Verb smalf enough, however, so that it coifld be assumed 
that the linear ranges were not exceeded. 

CONCLUDING REMARKS 

A flight investigation has been performed on a 35° swept- 
wing airplane in which the dynamic response characteristics 
were measured. Transient responses to elevator, rudder, 
and aileron deflections were recorded through the Mach 
number range of 0.50 to 1.04. The following remarks can 
be made regarding the longitudinal and lateral-directional 
characteristics. 



Moch number 


Fiuobb 16. — Variation of the derivatives Cng, Cn^ Ct^, and Ctg with 

Maoh number 



Mach number 

Figube 16. — ^Variation of the derivatives and with Maoh 

number. 
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tONGITTOINAL CHABACTBEISTICS 

Pitching velocity and nonnal-acceleration responses could 
be defined by simple second-order transfer functions as 
predicted from linear equations of motion. Except for the 
erratic variation in pitch damping between Mach numbers 
of 0.00 to 0.95, experimentally determined stability deriva- 
tives showed no unusual trends and agreed reasonably well 
with wind-timnel predictions. 

nATEHAL-DIEECTIONAl CHAHACTERISTICS 

Airplane responses in yawing velocity and sideslip angle 
duo to rudder disturbances were represented by second-order 
transfer functions that are related solely to the oscillatory 
mode. Simple first-order equations adequately defined the 
rolling velocity response to an aileron input. In either case 
the simplified equations closely defined the measured time 
histories while describing the frequency responses through 
the range of 1 to 10 radians per second. 

Fourth-order transfer functions calculated from the basic 
equations of motion using wind-tunnel and theoretical esti- 
mates of the various stability derivatives could he simplified 
by neglecting small quantities and by making approximate 
cancellations until they were of the same form as those 
evaluated from fiight data. Furthermore, it was possible to 
express the coefficients of these transfer functions in terms 
of individual stability derivatives. Frequency responses com- 
puted from these simplified equations were almost identical 
(between 1 and 10 radians per second) to those computed 
from the exact fourth-order transfer functions, and when 
compared with experimental resffits there was generally good 
agreement. Thus it is concluded that the simplified transfer 
functions form a reliable basis not only for estimating air- 
plane responses but also for the flight evaluation of stability 
derivatives, and the methods used here are felt to he suffi- 
ciently general to apply to any conventional airplaine. 

Experimental values of the derivatives C'npj 
(7„, , and (7ij compared favorably with predictions, baaed on 

theory and wind-tunnel measurements, at Mach numbers 
below 0.95, while at higher speeds, where predictions are 
questionable, there was some deviation. There were also 
notable discrepancies in flight values of <7?^ and <7#^, ob- 
tained at the 10,000 foot altitude, which were attributed to 
structural deformations resulting from aerodynamic and in- 
ertia loads. Additional evidence of aeroelasticity appeared 
when the Fourier analysis for was extended to fre- 
quencies beyond 10 radians per second; the frequency 
response showed evidence of aeroelastic deformation which 
appeared as an additional mode of motion not consistent 
with rigid airplane theory. 

Aitbs Abhonautical Laboratory ' 

National Advisory Committbb for AlBronautics 
Moffett Fibld, Calif., Sept. 17, 196S 
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APPENDIX A 

TRANSFER OF AXES 

The equations of motion normally used in airplane dy- 
namics are based on a system of axes fixed in the airplane 
m which the X axis is the intersection of the plane of sym- 
metry and a plane perpendicular to the plane of symmetry 
that contains the relative wind vector. These are normally 
referred to as stabihty axes. The angular displacement 
between the longitudinal stability axis and the reference 
axis of the airplane is equal to the trim angle of attack ar. 
Since recording instruments are generally alined with the 
reference axis, measurements of angular rates may he 
corrected to conform to stabihty axes notation as indicated 
in figure 17. Here p and r are vector components of the 
resultant rotation of the airplane, and the subscript 1 refers 



From figure 17 it can he seen that 

P=Pi cos ar-hri Sin. ar 
r=Ti cos ar — Pi sin ar 

Sideshp angles can he transformed by the relation 
’tan /3=s=tan /3i cos ar 

For most purposes these conversions are required only 
when the angle of attack is large. In this investigation it 
was found that the corrections to rolling velocity could he 
neglected in all cases because of the high roll to yaw ratio 
when using either a rudder or aileron input. In the cases 
of the yawing velocity records, however, the correction was 
sizsffile, particularly with tli^ aileron input. The correction 
to P was neglected in every case. 
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APPENDIX B 

CONTROL INPUTS 

Wlieii frequency responses are to be calculated from 
transient records, care should be given to the choice of a 
suitable forcing function. The frequency range through 
which accurate transformations can be obtained is definitely 
limited by the shape of the control input. Theoretically, 
a pure impulse (zero time dmation) is the most desirable 
input for all purposes because it gives uniform excitation 
to the entire frequency spectrum. The transform of a step 
input, on the other hand, has a magnitude that varies 
inversely with frequency and thus gives infinite excitation 
to the zero frequency component at the expense of the 
higher frequencies. 

The nearest physical approach to a pure impulse is an 
input that is roughly triangular in shape as shown in figure 18. 
Letting a equal the slope and T equal the time base of the 
triangle, the Fourier transformation of this input can be 
obtained from the relation 



Integration results in a transformation with the following 
real and imaginary parts: 


s-2f 

03 

osT / 
cos-^( 


(Bl) 

. 2a . 

osT / 

osT .\ 

(B2) 

sm 

_.COS— ly 


indeterminate points. This is shown in figure 18 where the 
transform magnitudes of two triangular pulses are plotted. 
One has T=1 second and a=4, while the other has T=% and 
tt=16. The areas under the two triangles are equal so that 
the transforms have equal magnitudes at zero frequency. Re- 
ducing T from 1 second to K second doubles the period and 
moves the first indeterminate point from a frequencj”^ of 
45 t radians per second to 8v. For pmposes of comparison, 
transformations of a imit step and a unit impulse are also 
shown. 

As T is further reduced, the magnitude of the triangular 
pulse more closely approaches the constant value that is 
characteristic of tiie pure impulse. To gain full advantage 
from the smaller T, the slope must be increased to maintain 
the same area under the pulse. A practical limitation is 
fixed by the maximum rate at which a control surface can 
be moved, and any further reduction in T results in smaller 
over-aU magnitudes. The most desirable input, therefore, 
is a compromise between large area and short time duration. 

From figure 18 it would appear that a pulse-type input is 
well suited for detennining low-frequency characteristics. 
However, the following explanation will show that this is 
not true. 

Generally it is impossible to return a control surface pre- 
cisely to its initial position after application of a pulse input. 
Even if a cbs.in stop or other device is used there is still apt 
to be a small residual deflection 5(T) after time T. 

If 5(T) is exactly zero, then as w approaches zero (from 
eqs. (Bl) and (B2)) 

J->0 



FiauBB 18. — ^Fourier transfona magnitudes of a step, an Impiilse, and 
two triangular-shaped inputs. 


The magnitude of the transformation is then 
III 2a/, coT\ 

It can be seen that 1 5 1 is periodic and is zero when u—ivIT, 
Sv/T, .... At these frequencies the transform of the re- 
sponse to this input would also be zero, and thus the ratio of 
output to input woidd be indeterminate. A reduction in 
T would increase the period and reduce the number of 




aT^ 


Now if 5 (T) is finite, the real and imaginary parts of the 
transform of the entire input are 

2a . wT { uT 5{T) rr 

I=—i sm ( cos 1 ) COS wi 

w* 2 V 2 / w 

„ 2a o)T ( ^ uT\ 8{T) rr 

R=—o cos ( 1— COS ) sin oii 

w* 2 \ 2 / oj 

In this case as oi— >0, 


/— >co 

R-^^-5(.T)T 

thus the zero frequency magnitude is infinitely large regai’d- 
less of how small 5(T) may be. Therefore, even though 8(T) 
appears to be zero on a flight record (i. e., 5(T) is less than the 
least count of the recording instrument), there is still the 
possibility of an infinite error at zero frequency. A step 
input is not subject to these large low-frequency errors; an 
error of 1 percent in the reading of the step deflection merely 
means an error of 1 percent in the transformation. 
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APPENDIX C 

PREDICTED AIRPLANE RESPONSES 


For constant forward velocity and only small disturbances 
from trimmed flight, the equations that define the motions 
of an airplane can be divided into two independent groups, 
the longitudinal equations and the lateral-directional equa- 
tions. For convenience, these are handled separately in 
developing theoretical transfer functions. 

LONGITUDINAl BESPONSES 

In view of the assumptions made, the longitudinal equa- 
tions may he written as 


mV(a-g) = -a^iCr,^a+ 

q_o8<^ 


(Cl) 


Dividing the equation (Cl) by mV and (C2) by J, the equa- 
tions may be rewritten as 


a — g=ZaCt-\-Zi^t 


(C3) 

(C4) 


Applying the Laplace transformation for za-o initial condi- 
tions and solving simultaneously for g and a gives 


where 


q C ieD-\-Oot 

S, D^+bD+k 

^ GiaD-\-Ooa 

h=-Za-M-M^ 

k=ZaM,—Ma 

Oot=MaZi^ — MsJZa 

a.=z>, 

Ooa=AISg — Zt^j 


(C5) 

(06) 

(07) 

(08) 
(09) 

(CIO) 

(011) 

( 012 ) 


The incremental change in normal acceleration is defined as 
Thus 


~ B^-VbD+k 


where 


Ckn=VZi, 

Con= - V{MaZ,-M>,Z.) = - VCa, 


(013) 

(014) 

(015) 

(016) 


In most conventional aircraft Ca„ and Om are very small 
compared to On and thus 

On 


n 


5f D^-\-bD-\~k 

The stabfiity derivatives Cn„, 0„^+Cn^, and may be 

obtained from approximate expressions for the transfer- 
fimction coefficients. 

The term can be determined from k (eq. (08)) by 
omitting the term ZaMq. Investigation has shown that, in 
tbipi case, tbia term is very small compared to Ma- Thus 


k=-Ma=-G„ 


qoSc 


or 


P =—J^h. 

qoSc 


(017) 


The dampiug factor 0m^-\- G„^ can be obtained from equa- 
tion (C7) by using known values of 6 and Ol^. 


M,+Ma=-Za-h 


or 


O -\-0 — 




(018) 


The elevator effectiveness derivative can be deter- 
mined from equation (09). The term is, for conven- 

tional aircraft, very small as compared with Ms,, and thus 


or 


Ms, — Oil, 


G — P 


(019) 


LATEBAL-DIBECnONAL RESPONSES 

Equations of motion. — ^The three equations that define 
the lateral and directional motions of the airplane may be 
written as: 

qoSb GtpP—qo8b Oi^d (C20) 

<p+{lJ)^-qoSb 

goSb 0„^=qoSb 0,,S (C21) 

(— TF cos i)<p-\-{mW—W sin y)f+ 

(mVD-q„SGpp=q,S0y,8 (022) 

By dividing equation (C20) by Ix, equation (021) by Iz, 
equation (C22) by mV, and by introducing new symbols, the 
three equations become: 

(IP-L,JD)v>+ (-rxD^-LrD)ik-L^P=L,8 (023) 
i-rzIP-N„D)v>+(T>^-NrD)^k-N^p=Ni8 (024) 
-Kiv+(.D-K,)>p+{D-Te)p=Yi8 (025) 
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CHABACTERISTIC EQUATION 

The characteristic equation A is formed by expanding the 
major determinant to give 

^=D{GJ>+OJJ^+CJ>^+OJ)+Oo) (C26) 

where 

C4=l — TxTz 

Cz— — Lp — Nt — Tfl(l — rrrz) — txN p — tJLt 
Gz= {Nr+TxLr) Fp+ (Lp+TxN;) (L^-NJ.r ) + 

(Nfi-\rrzLft) 

G,=- (L^r-NpLr) Fp+ {LfNp-N^,) -K, {Lf^+rxN^) - 

^3(Np+rzLp) 

Go= -KiiLrNn-NrLff) -KziLpNp-N^^) 

These coefficients can be further simplified by making the 
following substitutions: 

Let 

L/=Lp+rxNp N/=Nr+r^r 

Lff =Lfi-\-TxN N fi' =N ^-\- txL^ 


Then the coefficients of the characteristic equation are finally 
expressed as 

Gi=l — VxTx 

Gz=Lp'-N/-Y^CL-rxr^) 

Ga= (N/+L/) Y,+ (LpNr-NpLr) +N/ 

(7i= - (,LpNr--NpLr)Y»+ {L^p-N^p)-K,W-KzN/ 

a K^{L,N»-NrL^)-Kz{L^p-N^p) 

In factored form, equation (026) is 

A=D(jD — Xi) (P — Xj) (P — X3) (P — Xi) 

where Xi and Xj are designated as the spiral and rolling roots, 
respectively, and where X 3 and X< are a complex pair (v±iwi) 
that describe the oscillatory mode. 

For convenience in this investigation A has been expressed 
as 

A=D(P — Xi) (P — Xs) (Z?^-l-CiZ?-l-C2) 

where 

Ci= — (X 3 +XJ = — 2<r 
Ca = XsX^ = CO,* 


Here Ci and Cj are real coefficients that define the damping 
and period of the oscillatory motion. 

The quadratic term may also be written in the form 


where 

and 


\ C-J* V 


“b= vs 



SIMPIJFICA'nONS OP THE CHABACTEHISTIC EQUATION 


A=miy^- (L,+Nr+ Yi>)D^+ (LpY,+LpNr+Nf)D-LpNp] 

The cubic term can be factored exactly so that 

A=D^{D-Lp)[D^- (iV,+ Yp)D+Nf] 

This form of the characteristic equation considers only the 
oscillatory and rolling modes. It enables the coefficients 
Gt, Gi, and Cs to be expressed directly in terms of aerodynamic 
derivatives or simple combinations thereof. 

Even when the product of inertia is significant, the char- 
acteristic equation may be factored approximately into the 
simple form 

A=D^P-Lp)[D^- (Nr+ Yp)D+N0'] 

‘While the factorization is not exact, it is nevertheless 
justifiable in many cases. 

TBANSFEE FUNCTIONS 

From the three equations of motion the airplane responses 
in <p, yp, and j 8 can be readily calculated. In the following 
equations, 3 refers to either an aileron or a rudder disturbance. 

,p_AP^+AP^+AJD-\-Ao 

5 A 

where 

Ai=Li-\-rxNi 

Az= -UiX^+Nr) +Nz P-rxY^) + YPp' 

A,=U{NrY^+Np)-NiiU-^LrY^)+Yz{LrNp-NP^) 

Ao=KzQsfiL^-N^i) 

5 A 

where — 

Bz—Ni-\-rzLs 

Bz= -N,(Lp+ Fp) +LiiNp-rzYi,) 4- YiN/ 

Bi=NPpY^-L>NpYp+ Yi(LpNp-NfiLp) 

Bo=Ki(LzNp-NP^) 

PjDiEzD^+EjD^+EiD+Ep) 

d A 

where 


Et=Yi(l—rxrz) 

Ez= - Yz{lY/+Lp')-Lirz-Nz 

E^=YzP^-NpLr)-U{Np-rzKz-K{)+Ni{rxK,+ 

Lp-^K^ 

Eo=K^{NP-LzNr)+Kz(PJSfp-NPp) 


Simplifications can also be made in these expressions by 
neglecting small quantities; however, this can be shown more 
clearly in the numerical example that follows. 


If Go and other small terms are neglected when the flight- 
path angle y and when the product of inertia are essentially 
zero, the characteristic equation may be written as 


NUMEEICAI, EXAMPLE AT A/- 0.8 

By use of values of stability derivatives shown in table H, 
and with 7 = 0 , predicted responses for p=D(p,r=DTp, and ^ 
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for both aileron and rudder inputs have been calculated and 
found to be 

A=D(IH+3.652l7+15.16l»»+41.26l>+0.0289) 

=Z?(Z?+0.00070) (I>+3.078) (Z>»+0.573l?+ 13.40) 

r_-7.60I>(D+3.091)(I>*+0.0270l>+0.208) 
dr A 

p_6.16Z)*(I>+4.436)(Z?-6.210) 

Sr A 

/3_0.0339D(Z>+3.053)(I>+225.3)(Z?-0.00703) 

Sr A 

r_O.099D(D+3.978)(I>“-1.758I>+7.358) 

Sa^ A 


p_Zp'D+Lr7Vff-NrZg 
^ Na'D-LrJSff>+Nr,L^ 


The quantity LrNa — NrL^ in the numerator is generally 
very small and can be neglected. In the denominator, the 
term is small compared to the other terms and can also 


be neglected. 


The ratio, is then simplified to 


<p La'D 

T~Na'{D-L,) 


When the complex root Xs=(r+'i«i is substituted for the 
operator D, this expression is the ratio of the free oscillatory 
responses of </> and V' at any time t. The ratio of p to r is 
obviously the same and can be expressed as 


p_36.4D*(I>»+0.655Z?+13.68) 

5a"" A 

/3_0.0008D(D+0.990)(Z)-1.094)(D-870) 

5a^ A 

When small terms are neglected, r/5, can he expressed as 


p_ La'((T+ica{) 
r Nfi 

The actual magnitude of this ratio is 

p La 


r_ -7.60Z)2(D +3.091) 

Sr D*(I>+3.078)(I>*+0.573D+13.40) 

and then by an approximate cancellation this reduces to 

r_ -7.602? 

5, 2?»+0.6732?+13.40 

Similarly, fi/Sr can be simplified by neglecting small terms so 
that 

/3_ 7.642?’(2?+3.053) 7.64 

5r""2?*(Z?+3.078)(2?»+0.6732?+13.40) ”2?>+0.573Z?+13.40 

It can be seen that this expression for fi/Sr is practically 
identical (with opposite sign) to the integral of the simplified 
equation for r/5,. It is also possible to simplify p/Sa as 
follows: 

p_ 36.4Z?*(2?*+0.6552?+13.68) 36.4 

5a""2?*(2?+3.078)(2?»+0.6732?+13.40) ^2?+3.078 

Similar simplifications have been made for other Mach 
numbers and found to be equally valid. 

APPENDIX D 

FLIGHT EVALUATION OF Cia 

From time histories of the free oscillatory responses of p 
and r, it is possible to evaluate the derivative Ci^ as follows. 
If the rolling moment and yawing moment equations (C23) 

and (C24) are set equal to zero, then the ratio ^ of the free 

oscillation may be obtained by simultaneous solution of the 
two equations to give 

¥)_(rx2?»+ LrD)Na+(D^-NrD)Za 
^ (D^-t^)Na+(rzD^+N,D)U 


In this form \p/r\ is the ratio at any time t of the amplitudes 
of the envelopes that enclose the oscillatory motions of p 
and r; <r is the rate of damping of the envelope; and wi is the 
natural frequency of oscillation. When <r is very small as 
compared to it can be omitted; thus 


La'ui 

La' 




If Na' and Lp are known, it is then possible to evaluate 
La from measured time histories of p and r. 
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TABLE I.— PHYSICAL CHARACTERISTICS OF TEST AIRPLANE 


Wing 

Total area 287.9 sq ft 

Span 37.1 ft 

Aspect ratio 4.79 

Taper ratio 0.51 

Mean aerodynamic chord 97.03 in. 

Dihedral 3“ 

Sweepback of quarter-chord line 36° 14' 

Aerodynamic and geometric twist 2° 


Root airfoil section (normal to quarter-chord NACA 0012-64 (modified) 
line). 

Tip airfoil section (normal to quartenehord line) _ NACA 0011—64 (modified) 
Ailerons 


Area, each 18.6 sq ft 

Span 9.18 ft 

Chord, average 2.03 ft 

Defiection, maximum 14° up, 14° down 

Inboard end at 61.6 percent 6/2 

Horixontal tail 

Total area (including 1.20 sq ft covered by 36.0 sq ft 
fuselage). 

Span 12.8 ft 

Aspect ratio 4.66 

Taper ratio 0.46 

Main aerodynamic chord 34.7 in. 

Sweepback of 0.25 chord line 34°36' 

Airfoil section NACA 0010—64 

Tail length 18.12 ft 

Vertical tall 

Area, total 34.4 sq ft 

Span 7.6 ft 

Aspect ratio 1.74 

Taper ratio 0.36 

Sweepback of quarter-chord line 35°00' 

Rudder 

Area 8.1 sq ft 

Span 6.6 ft 

Chord, average 1.23 ft 

Defiection, Tnavimiim 24.8° right, 26° left 

Take-off weight (o. g. at 23.0 percent MAC) 14,102 lb 

Tending weight (c. g. at 21.9 percent MAC) 11,612 lb 

Average weight for calculations (c. g. at 22.6 percent 12,800 lb 
MAQ. 

Moment of inertia about X axis 7,246 slug ft’ 

Moment of inertia about Y axis 17,480 slug ft’ 

Moment of inertia about Z axis 23,190 slug ff^ 

Inclination of principal longitudinal axis with respect — 2.6° 
to fuselage reference axis. 

















